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Boots of a "Folynooial (DOI oi-ly) 

Bitlre program 

Depends on condition of the polynomial. Uaually about 9 

decimal places . 

0.082 n + 5.2 n + 18 seconds 

Tina time will depend partly on the distribution of roots 

This program calculates tlie roots of a polynomial 



mt \ * a-1 . 

f (x) • m Q x + a x x + — 



— + a. 



vmose coefficients mavo boea patched on tape. Am m'tb 

and n complex root* obtained. Real numbers are treated as 
complex numbers vhose imaginary parts are aero. 
The master tape is read into the computer and then tho tape 
containing the coefficienta is read. As each root is com- 
puted it will be punched . After all n roots have been punched 
the computer will stop on a black switch stop*and another 
tape of coefficients may be read. Hence one may prepare 
several equations on a single tape and have them all solved 
by leaving the black switch in the disable position. If this 
is done the tape should be terminated with the symbol N. 
Whoa this symbol is read in place of a new set of coefficients 
the routine will stop on OF. 
PREPARING TIE TAPE OF COEFFICIENTS To prepare a tape of coefficients of 

the polynomial 



n n-1 







ene should first type n as a deciral integer followed by a 

carriage return. The n + 1 coefficients are now typed in 



order a 



o' V 



— -. a 



n 



Each must be followed by a carriage 
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return. Each coefficient a Must Toe regarded as a complex number in floating 

decimal form 

a . (A + iB) x 10* 

where |A| < 1 and |b| < 1 and p is an integer. 

The three nuatbera A, B, and p are typed in order, A and B as 
signed fractions vith the decinal point preceding the first digit in each case, 
and p as a signed integer. The absence of digits after a sign is treated as 
zero. 
EXAMPLE The equation 

x 5 + 2kx k + (3 - 6Hi) x 3 - (.05 + .003*1) x 2 ♦ .39 - 

would be prepared as 

5 

+1++1 

+2W+2 

+03-#H-2 
.5-03^-1 

+++ 

+.39*+ 
FORM OF THE MSULTS The roots are printed in floating decimal fonA B p, 
with A and B printed as signed fractions to nine decimal placet and p printed 
as a signed integer. Each root Is followed by a residual printed to 3 decinal 
places which is obtained by substituting the root back into the original equation. 
Tke residuals also appear In -floating decimal form. 

MATHEMATICAL METHOD The program u?ea a 1.8th order iterative procedure to 
find a root. Each successive iterant is obtained by solving for the nearer root 
of the quadratic passing through the last three iterants . This solution is 
accomplished by a variation of the standard quadratic f arasula . By use of this 
method the dumber of evaluations of the f- \ction is reduced, and hence the time 
to solve equations of high degree. As eai.'.i /not is found it is divided into the 
polynomial, thus reducing its degree by one, and tfc process is repeated. 
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as follows 
forault 



The specific algorithm which is used my toe described 
Given quantities f Q , f., fp, x, X, 5, h we calculate V by the 



-2f 2 5 






5 X (f n X - t . b + f ) 



where b * f Q X 2 - f . o 2 + f (X + 6) 

The sign before the radical is chosen so as to Wkkm 
denoednator have the larger Magnitude. In one iteration we replaca 

X by X' 
h by X'h 
xbyx+ X'h 
5 by 1 + X 1 
f 2 by f (x + X ! h) 

f i * f a 

1 by t x 

Before calculating each root ve start with the initial 



values - 



f Q - f(0) - f f (0) + 1/2 f"(0) 
t l « f(0) + f'(0) + 1/2 f"(0) 

t 2 = f (o) 

x - 

X m - 1/2 

6 =» + 1/2 

h « - 1 

Ve make use of the relations 

f(0) . a n , f»(0) = a n ^, 1/2 f"(0) - a n _ 2 

A final value of x, say x., is printed when 



i x i - x i-i 



.9 
< 10 , that is when the change in x is 



no more than one part in 10 



k - 



In order to ensure convergence a special process is used 
whenever f /f > 10. ;n this case \' iz 'replaced by l/? V and h, x, and f" 2 
are recomputed . The original process is then resumed. Final convergence is 
not- affected by this process . 
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Convergence criteria* 
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